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I  -  LOW  TEMPERATURE  QUARTZ  CRYSTAL  OSCILLATOR 


GENERAL  PRESENTATION 


1/f  noise  remains  the  main  limitation  to  the  stability  of  quartz  oscillators, 
when  thermal  noise  and  random  noise  are  reduced.  It  has  been  demonstrated  that 
1/f  noise  level  depends  on  the  resonator  Q-factor  following  a  1/Q4  law. 

Q-factor  has  a  finite  value  on  account  of  phonon  interactions,  but  also  because 
of  impurity  relaxation  and  diffusion.  Therefore  in  this  second  part  an  analysis 
of  impurity  mechanisms  was  undertaken  to  determine  their  contribution  to  Q-factor 
and  frequency  fluctuations  at  low  temperature. 

Preliminary  results  on  short  and  long  term  frequency  stability  of  quartz  reso¬ 
nators  in  liquid  helium  were  obtained  by  using  an  adapted  dual  crystal  passive 
system.  The  electronics  was  operated  at  room  temperature  and  the  crystal  was 
in  a  copper  can.  Temperature  was  first  controlled  at  4K,  and  then  near  1.5K  in 
superfluid  helium. 


LOSS  MECHANISMS  DUE  TO  IMPURITIES 


1 )  Relaxation  of  interstitial  alkali  ions 


Fig.  1 

Double-well  potential 


The  most  common  interpretation  of  ionic  con¬ 
ductivity  in  quartz  crystals  assumes  that  some 

4+  3+ 

Si  atoms  are  replaced  by  A1  associated 

with  a  monovalent  alkali  ion  (Na+,  Li+  or 

K+  ..)  achieving  change  neutrality  in  the 

( 2 ) 

crystal  .  The  ion  occupies  an  interstitial 
site  and  is  trapped  in  a  double  well  poten¬ 
tial  (Fig.  1).  It  can  jump  from  one  position 
to  the  other  one  by  thermal  activation  and 
cross  over  the  potential  barrier  V  . 


The  impurity  is  characterized  by  two  quantum  states  with  energy  levels  E  , 
and  energy  splitting  E  =  E  -  E  .  If  and  are  the  probability  densities 
of  transition  respectively  from  state  1  to  2  and  2  to  1,  in  the  absence  of  sound 
wave,  a  relaxation  time  x  can  be  defined  as 


1 

x 
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(1) 


If  a  sound  wave  is  superimposed  the  energy  spliting  E  is  modulated  by  the  cor¬ 
responding  high  frequency  strain  n  and  takes  the  instantaneous  value 

E'  =  E  +  Dn  (2) 

where  D  is  the  deformation  potential. 

Let  nQ  be  the  number  of  impurities  per  unit  volume,  and  the  probability  of 
the  impurity  to  be  in  state  1  with  energy  E^  when  in  equilibrium. 

Px  =  ( e8E  +  l)"1  (3) 

with  8  =  1/k  T.  T  is  the  absolute  temperature  and  k  the  Boltzmann  constant. 

D  D 

(  3  ) 

The  energy  contribution  of  the  impurities  can  be  written 

Et  =  no  ^P1  (E  +  Dn)  +  E2^  (4) 

Differentiating  with  respect  to  strain  gives  the  impurity  contribution  to 
the  elastic  stress. 


(5) 


n  D*n  <5P 

.  _  _o _ 1_ 

CT  ~  1-iuT  6 ( E ' ) 

where  w  is  the  angular  frequency  of  the  sound  wave. 


This  relation  shows  that  when  the  impurity  system  is  disturbed  by  the  sound 
wave,  a  relaxation  process  occurs  with  time  constant  t  to  restore  equilibrium. 
There  is  a  feedback  of  the  change  in  population  of  impurity  states  to  the  elas¬ 
tic  stress  which  leads  for  the  sound  wave  to  dispersion  (feedback  in  phase)  and 
to  absorption  (feedback  out  of  phase). 


The  ratio  <5c/n  gives  the  change  <5C  of  the  elastic  constant,  from  where  is  easi¬ 
ly  derived  the  velocity  change  5v  and  the  attenuation  a  of  the  sound  wave 


<5V  =  — —  &  (SC) 
2pv  e 
o 


a  =  -  — r  I  (« C) 

2pv  m 
o 


(6) 

(7) 


p  is  the  specific  mass  of  the  crystal  and  vq  the  sound  wave  phase  velocity. 

From  these  two  relations  can  be  obtained  the  frequency  shift  and  the  Q-factor 
of  the  resonator 


,  n  D  (i)  t 

1  _  O _  ,  o  , 

0  4  D  V*  k  T  4  +  w2  T*; 

o  B  o 


OJ  —  0) 


n  D2 


(t— TT-rr) 


8  p  v2  k  T  '1  +  w2  t  : 
o  B  o 


(8) 


(9) 


The  process  which  allows  the  impurity  to  cross  over  the  potential  barrier  Vq 
of  the  double  well  potential  is  a  thermal  activation  process.  It  can  be  easily 
shown  that  it  is  governed  by  an  Arrhenius  type  law 
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2v 
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-  V  /< E  > 
o  o 
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( 10) 


where  vq  is  the  jump  frequency  of  the  impurity,  and  <Eq>  the  mean  energy  of 
the  harmonic  oscillator  composed  of  the  impurity  in  one  of  the  two  quadratic 
potential  wells. 


k 


which  does  not  depend  on  the  deformation  potential. 


By  substracting  in  the  curves  of  Fig.  2  the  background  losses  due  to  phonon  in¬ 
teractions  the  contributions  of  impurity  relaxation  were  obtained  and  are  given 
on  Fig.  3. 


Adjusting  parameters  of  relation  8  theoretical  curves  (dotted  lines)  were  fitted 
to  the  experimental  ones  (continuous  lines).  This  curve  fitting  enables  one  to 


evaluate  the  value  of  the  relaxation  time  constant  t 
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5.10  s,  the  potential 


barrier  height  Vq  =  0.040eV  and  the  distance  between  the  two  wells  2d  =  8  A. 


These  results  show  that  the  low  temperature  relaxation  gives  information  or  he 
Al^+-M+  defect  structure.  The  impurity  content  nQ  cannot  be  directly  obta:  d 
from  such  measurements  because  the  losses  1/0  depend  on  n^D2  where  the  del  „a- 
tion  potential  D  is  not  known.  However  n^  could  be  determined  by  comparis 

with  a  reference  crystal  if  its  impurity  concentration  were  known  and  if  th. 
cut  and  the  vibration  mode  are  the  same. 


If  one  considers  that  the  relaxation  time  constant  t  is  fluctuating  by  6t,  this 
will  produce  a  fluctuation  6w  of  the  resonance  frequency.  From  relation  (9) 


6(ij 

ui 

o 


AE 

E 
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with  E  =  pv2 
o  o 


( 1  +  (i) 4 


and 


AE  = 


_  _6r_ 

>*  T 

n  D2 
o 

4kBT 


(14) 


Using  relation  (8)  in  (14)  gives 


This  will  also  exhibit  a  1/Q  -law  between  frequency  and  relaxation  time  fluc¬ 
tuations  in  the  power  spectral  densities.  This  was  already  pointed  out,  with 
the  same  form,  in  the  case  of  fluctuations  due  to  phonon  interactions  ' ^ . 


2 )  Sound  wave  phonon  scattering  by  static  defects 

In  addition  to  the  previous  relaxation  of  interstitial  alkali  ions  is  the  scat¬ 
tering  of  the  sound  wave  by  the  substitutional  aluminium  atoms  replacing  sili¬ 
con.  This  can  be  considered  as  an  isotropic  effect  which  induces  alteration 
of  the  kinetic  energy  by  mass  difference  and  of  the  potential  energy  by  change 
of  the  linkages  (equivalent  to  a  change  of  some  elastic  constants).  Both  lead 
to  similar  absorption  and  only  scattering  by  atoms  of  different  mass  will  be 
considered  here. 


A  linear  chain  of  atoms  with  mass  m  has  some  substitutional  atoms  of  mass 
m+Am  ;  the  perturbed  Hamiltonian  can  be  written 


W  =  W  +  AW 
o 


where  Wq  corresponds  to  a  perfect  crystal  and 


AW  =  V2  Am  (3u  /3t): 

n 


u  is  the  displacement  of  an  atom  at  the  lattice  site  R 
n  n 


,  jk.R  -jk.R 

1  r  1/  ti  r  n  +  n. 


(16) 


(17) 


(18) 


The  sum  is  taken  over  all  normal  modes  with  wave  number  k  and  frequency 
G  is  the  number  of  atoms  in  the  crystal,  a^  and  a *  respectively  are  the  annihi¬ 
lation  and  creation  operators  of  the  harmonic  oscillator.  Using  (18)  and  (17) 
it  can  be  shown  that  for  n^  impurities  per  unit  volume  randomly  distributed  in 
the  crystal 

n  Am  ti 
...  o 


l\f 

kk' 


0)  U) 

k  k' 


+  a, 


( 19) 


Let  et  N  be  the  phonon  numbers  in  normal  mod's  k  and  k'.  Applying 
Fermi's  Golden  Rule  gives  the  transition  probability  and  then  the  rate  at  which 
phonons  are  scattered  from  mode  k. 


dN, 

_ _ k 

dt 


n  " 
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Am' 


k  k' 


'<V1)Nk' 


-  VV*1 


.  U)  —CO 

'  k  k' 


-  7  - 


(20) 


6  is  the  Kronecker  symbol  and  6 (w.  , -io  )  assumes  conservation  of  energy.  If  con 

K  K 

sidering  that  all  modes  k1  are  in  thermal  equilibrium,  equation  (20)  can  be 
easily  solved  However  if  interactions  between  thermal  phonons  are  taken 

into  account  a  system  of  first  order  linear  dil  '.‘rential  equations  must  be  con 
sidered  as  a  whole. 


Let  n  and  n  be  the  phonon  numbers  when  in  thermal  equilibrium 

-  ir‘ 
nk  -  ie  -  11 


\  =  "k  *  nk 

V  ■  V  *  V 


The  conservation  of  energy  implies  w  =  w  and  therefore  n 

K  K  K 

set  of  equations  similar  to  (20)  becomes 


n,  ,  .  Thus  the 
k’ 


K.  ! ,  X  , 

Hr  "  K  L  (nir<  ~  nL 


n  if  .  z 

with  K  =  (— )  u) 2  and  w 

8G'  m  o  o 


w.  =  to 
k  k' 


Time  dependent  solutions  of  the  form  =  a^  e  can  be  used.  This  leads  to 
the  new  time-independent  system 

ak  *  K  T  l  (ak  ~  akJ  (25) 

k' 

which  has  the  eigenvalues 


f—  -  (p  +  i: 


p  being  the  total  number  of  possible  values  of  k'  at  the  frequency  The  de¬ 
generated  eigensolution  is  equivalent  for  large  p  to 

1/Kt  »  p  (28) 


The  total  number  of  k'  values  is  obtained  by  integrating  in  the  k-space  over 

the  surface  w  =  u> ,  ,  =  w  and  taking  into  account  the  three  acoustic  branches 
k  k '  o 


8 


, V . V o  .*■  .v yv 


Ls  u>* 
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where  3/C3  =  2,/C3  +■  1/CS .  C  and  C  respectively  are  the  transverse  and  longi- 

o  u  %  X, 

tudinal  sound  wave  velocities  and  Ls  is  the  crystal  volume. 


3  n  L3 
o 


16it  C3 


Introducing  the  lattice  parameter  a,  the  resonator  0-factor  is 


3  n  a' 
o 

16tt  G 


This  shows  that  1/0  is  temperature  independant  and  will  give  a  threshold  for 
the  lowest  temperature  when  phonon  interactions  and  interstitial  impurity  rela 
xation  are  minimized.  But  this  effect  appears  to  be  neglegible  at  low  frequen- 
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LOW  TEMPERATURE  DUAL  CRYSTAL  PASSIVE  SYSTEM 


To  perform  frequency  stability  measurements  on  a  resonator  at  low  temperature 

with  the  electronics  at  room  temperature  the  most  convenient  system  is  the  dual 
(7) 

crystal  system  because  of  the  distance  between  the  inner  part  of  the  dewar 
and  the  outside. 

1)  Description 

In  fig..  4  is  shown  a  schematic  diagram  of  the  system  adapted  to  low  temperature 
measurements. 


Fig.  4 

Schematic  diagram  of  dual  crystal  passive  system 


The  output  signal  of  the  auxiliary  quartz  oscillator  is  phase  modulated  at  a 
fixed  frequency  and  then  is  passed  through  a  directional  coupler  to  the  quartz 
resonator  under  test. 

The  reflected  signal  by  the  quartz  is  amplified  and  then  detected.  The  auxilia¬ 
ry  quartz  oscillator  is  locked  on  the  resonance  frequency  of  the  quartz  inside 
the  dewar  with  the  help  of  a  lock-in  amplifier  which  is  operated  at  the  modula¬ 
ting  frequency  f  . 


The  heart  of  this  apparatus  is  the  phase  modulator  which  must  be  designed  with 
care.  As  a  matter  of  fact,  the  spurious  amplitude  modulation  generated  at  the 
same  time  as  the  phase  modulation,  must  be  minimized,  because  this  phase  modu¬ 
lation  after  reflection  on  the  quartz  gives  the  amplitude  modulation  the  modu¬ 
lation  index  of  which  is  proportional  to  the  difference  frequency  between  the 
auxiliary  oscillator  and  the  quartz.  The  set  constituted  by  the  directional  cou 
pier  and  the  reference  quartz  resonator  works  like  a  discriminator. 

_5 

It  appears  that  a  modulation  index  of  10  (which  gives  an  amplitude  of  -100  dB 

-13 

with  regard  to  the  carrier)  produces  a  relative  frequency  shift  of  10  .  The 

phase  modulation  is  an  original  design  and  is  entirely  digital. 


The  phase  modulator 

The  phase  modulation  is  obtained  by  using  the 
time  delay  produced  by  a  logic  gate  (Fig.  5) 
between  the  output  and  the  input  signals. 

Let  AT  be  this  time  delay.  It  corresponds  to 
a  phase  shift  A<t>  =  2*  AT/T  where  T  is  the 
period  of  the  HF  signal.  By  using  n  gates, 
it  is  possible  to  obtain  a  phase  shift 
between  0  to  nA$  with  n-steps. 


AT  T 
Fig.  5 


A  four  bit  multiplexer  (Fig.  6)  using 
eight  gates  leads  to  an  output  signal 
having  a  triangular  phase  modulation 
with  a  peak  to  peak  phase  shift  of 
2n/3  (Fig.  7). 

This  modulation  waveform  is  obtained  by 
connecting  together  the  multiplexer  inputs 
1-15,  2-14,  3-13,  4-12,  5-11,  6-10,  7-9. 
The  fundamental  frequency  modulation  is 
equal  to  1/16  of  the  driving  frequency 
of  the  multiplexer  states.  This  driving 
frequency  is  obtained  by  dividing  the 
output  signal  of  the  auxiliary  locked 
oscillator. 
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Fig.  6 
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In  Fig.  3  are  respectively  represented  the  input  signal  (a),  the  output  signal 
of  the  first  gate  with  a  pnase  shift  of  AT  (b),  the  output  signal  of  the  second 
gate  with  a  shift  of  2AT"  (c)  the  command  signal  (d)  and  the  phase  modulated 

output  signal  (e). 


a)  5  MHz  input 


b)  Output  1st  gate 

c)  Output  2nd  gate 

d)  command  signal 


e)  phase  modulated 
output  signal 
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I  Fig.  8 

The  last  pulse  width  at  the  commutation  time  is  not  equal  to  Tq 

where  Tq  is  the  V2  period  of  the  5  MHz  signal.  In  fact,  when  the  phase  jump  is 

i  positive  (for  the  0  to  8  multiplexer  states)  the  pulse  width  increases  by  AT  and 

when  the  phase  jump  is  negative  (for  the  8  to  15  multiplexer  states)  the  pulse 

width  decreases  by  AT.  In  Fig.  9,  we 

show  the  pulse  width  versus  time,  we 

I  obtain  eight  pulses  larger  and  eight 

pulses  smaller  than  the  mean  T  .  Since 

o 

the  pulse  amplitude  is  constant  the 
index  modulation  can  be  expressed  by 


T+Ar 

o 


T-AT 

O 


pulsewidth 


i  m  1 1 1 1 1 


m  1 1 1 1 


time 


_ _ 8  AT _ 

Tqx  number  of  pulse  on  a  half  period  of  modulation 


Fig.  9 
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For  a  400  Hz  frequency  modulation  with  a  AT  of  8  ns  and  a  fundamental  frequency 
of  5  MHz,  this  modulation  "index"  is  equal  to  6  x  10-5  (or  -83  dB). 

In  fact  experiments  give  -60  dB  ;  this  difference  is  due  to  the  overshoot  when 
the  commutation  occurs. 

The  amplitude  modulation  rejection  factor  can  be  increased  by  using  a  monosta¬ 
ble  which  gives  a  constant  pulse  width. Thus  the  amplitude  modulation  can  be  re¬ 
jected  to  -100  dB. 


Buffer 

This  "cascade"  buffer  gives  the  possibility  of  obtaining  three  outputs  from 
one  input  with  a  feedback  coupling  of  100  dB.  Each  output  feeds  respectively 
the  phase  modulator,  the  divider,  which  gives  the  modulation  frequency  and 
the  third  output  which  gives  the  5  MHz  signal. 


Directional  coupler 


The  coupling  factor  of  the  directional  coupler  is  10  dB.  This  is  a  commercial 
one.  Its  directivity  is  40  dB.  A  better  directivity  would  be  desirable  but  is 
not  available. 


Amplifier  -  Detector 


It  permits  on  the  one  hand  amplification  of  the  signal  reflected  by  the  quartz 
to  a  level  sufficient  to  be  detected  by  a  diode  and  on  the  other  hand  presenta¬ 
tion  of  a  constant  load  to  the  output  of  the  directional  coupler,  which  is  the 
load  as  seen  by  the  resonator  under  test. 


Lock-in  amplifier 


This  lock-in  amplifier  consists  of  two  amplifiers  with  respective  phases  which 
are  out-of-phase.  The  outputs  of  the  amplifiers  are  sampled  at  the  half  period 
of  the  modulation  frequency  and  then  summed  on  a  resistor. 

The  sampling  signal  is  always  in  quadrature  with  respect  to  the  excitation  si¬ 
gnal  for  all  the  modulation  frequencies. 


Synthesizer 


The  resonance  frequency  of  the  quartz  at  low  temperature  presents  a  shift  of 
5  kHz  with  respect  to  its  value  at  room  temperature.  Therefore  the  auxiliary 
oscillator  cannot  be  a  Xtal  oscillator,  and  a  synthesizer  must  be  used  in  the 
loop.  This  enables  one  to  obtain  a  5  MHz  output  frequency  whatever  the 
frequency  of  the  resonator  is  in  a  4  to  6  MHz  range. 


Frequency  measurements  system 


The  first  FDM  was  built  by  the  LPMO.  Selective  filters  are  inserted  at  the 
output  of  the  frequency  multipliers  (x9  reference  channel  and  xlO  oscillator 
under  test  channel) in  order  to  prevent  the  beat  of  the  side  bands  of  the  two 
channel  signals  in  the  mixer.  Just  a  selective  filter  after  the  mixer  does 
not  avoid  these  spurious  signals,  because  they  are  at  a  frequency  near  the 
5  MHz  ±  10  af  desired  frequency. 

The  second  FDM  is  a  commercial  unit  which  enhances  by  a  factor  of  2x10^  the 
the  difference  frequency  of  the  two  signals  connected  at  its  inputs. 

The  measurement  system  noise  versus  integration  time  t  is  given  by  fig.  11. 


-  Frequency  stability 


The  frequency-temperature  characteristics  of  quartz  resonato  ah  low  tern 
perature  were  measured  between  4  K  and  1.5  K  and  are  given,  for  two  resonators, 
in  Fig.  13. 

AT-cuts  exhibits  a  general  behavior  corresponding  to  the  lower  curve,  where  the 
frequency  temperature  coefficient  goes  to  zero  with  temperature. 


Fig.  13 

Frequency  vs  temperature  for  FC  and  AT  cut 


Temoerature  (°K) 


This  is  wiidt  can  be  expected  from  the  theory.  At  1.5  k  the  temperature  co- 

_g 

efficient  is  of  the  order  of  6.10  /°K.  This  gives  an  indication  of  the 

-  5 

temperature  stability  which  is  necessary;  for  instance  1.6.10  K  to  achieve 

-11 

frequency  stability  of  1.10  . .  These  numbers  are  to  be  considered  only  for 

the  static  thermal  behavior  of  the  crystal.  The  dynamic  temperature  coeffi¬ 
cient  of  quartz  resonators  has  not  been  evaluated,  but  would  need  to  be  con¬ 
sidered  in  addition. 

In  the  present  state  of  the  research,  short  term  and  long  term  frequency 
stability  measurements  were  not  yet  performed  with  temperature  controls 
operated  at  the  same  time. 


Preliminary  results  were  obtained  with  a  resonator  in  the  helium  bath,  but 
without  temperature  and  pressure  control.  In  this  condition  the  result  in¬ 
dicates  a  frequency  stability  of  4.10-12  for  10  <  t  <  100  s  at  4  K. 

These  frequency  fluctuations  can  be  related  to  the  temperature  fluctuations 

-4 

of  the  free  running  bath  which  correspond  to  2.10  K  and  at  4  K  the  frequency 

—8 

temperature  coefficient  of  the  crystal  is  2.10  /°K. 


CONCLUSION 


Impurities  can  have  two  effects  at  low  temperature  :  diffraction  due  to  substi 
tutional  aluminum  atoms  and  relaxation  of  interstitial  alkalin  ions.  If  the 
last  one  seems  to  have  negligible  influence,  at  least  at  low  frequency,  the 
first  one  can  induce  important  relaxation  peaks  and  lead  to  frequency  fluctua- 

4 

tions  following  the  1/0  law  previously  pointed  out  in  the  case  of  phonon- 
phonon  interactions.  Experimental  verification  will  be  the  task  of  the  next 
part  of  the  project. 


Electronics  for  operating  quartz  resonators  at  low  temperature  and  for  control 
ling  temperature  fluctuations  were  developed  during  the  present  p.eriod  of  the 
project . 


-5 

Temperature  stability  of  2.5  10  K  was  obtained  at  1.5  K.  Taking  into  account 

_9 

the  frequency  temperature  coefficient  of  the  crystal,  6  10  K,  this  would  cor- 

-13 

respond  to  a  frequency  stability  of  1.5  10  .  Measuring  short  and  long  term 

frequency  stabilities  will  be  the  purpose  of  the  next  part  of  this  study. 
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II  -  FAST  WARM-UP  SAW  OSCILLATOR 


01  M  KM  I’KI  S|-‘\'i  AT  ION! 


Generally,  thermo i  behavior  of  SAW  oscillators  is  due  at  the  same  t  in"  t •>  : ; i ow 
temperature  variations  ( temperature  drifts)  and  fast  temperature  var  i  a  ‘  i  ,ris 
l  temperature  instabilities  of  thermostats  or  th'»rma  1  transient.:,  of  v,n:  -up 
r i 1 lators ) . 

!n  the  present  wort ,  eniy  fast  temperature  variations  are  t  o  be  con: . i  dep«.  q .  in 
this  rase,  a  spatial  distribution  of  temperature  within  the  crystal  arise.;. 
Thermal  effects  involve  temperature  gradients  which  induce  stress  and  strain 
distribution.  This  is  the  "Dynamic  thermal  behavior"  Thermal  stresses 

and  strains  are  coupled  with  the  high  frequency  surface  acoustic  wave  by  crys¬ 
tal  nonlinearities.  Therefore  velocity  and  frequency  are  shifted. 

In  a  previous  interim  report  (July  1981),  a  one-dimensional  model  of  temperature 
distribution  was  presented.  The  perturbation  method  used  to  calculate  correspon¬ 
ding  frequency  shifts  led  to  numerical  results.  Comparison  with  experimental 
data  showed  some  differences.  Consequently  a  new  model  of  temperature  distribu¬ 
tion  is  required  in  order  to  explain  the  dynamic  thermal  behavior. 

The  present  work  shows  that  : 

% 

-  the  measured  dynamic  temperature  coefficient  a  is  small  for  quart./  eats  having 
■ero  first  order  static  temperature  coefficients. 

-v  -(') 

-  experimental  values  nf  a  for  ST-cut:  and  l-ST-cut.  are  ‘,.f,  ,h  s/K  and 

1  !  .n-t.i'Hl'  ^  s/K  respectively. 

-  theoretical  results  obtained  by  using  a  new  one-dimensional  thermal  model  are 
in  good  agreement  with  the  experimental  results . 


The  theoretical  study  is  developed  along  the  following  lines  : 

-  calculation  of  temperature  distribution  in  a  om— d  i  men:;  iona  1  mode!  Inking 
into  account  fast  thermal  exchanges  between  the  crystal  plate  and  tie-  sur¬ 
rounding  medium. 

-  calculation  of  corresponding  stresses  and  strains  leading  to  the  t  h'-or'-tie.i 

'V, 

values  of  dynamic  temperature  coefficient  a  by  us  i  ng  a  per  t.urha  *  .<  on  met  hod. 

-  to  take  into  account  the  mechanical  mounting  of  the  quart  z-p  1  a !  e ,  .1  iwo-di- 
mensional  model  of  temperature  diffusion  is  developed. 


ONE-DIMENSIONAL  MODEL 


A  one-dimensional  model  has  been  used.  It  consists  of  a  crystal  plate  of  2h- 
thickness  and  of  infinite  length  along  Oy  and  Oz. 

x 


F  i  g .  1 

One-dimensional  tempera  turf*  d  i  s  tri  hid  i  >n 

The  bottom  of  the  plate  is  submitted  to  temperature  variation  $ft).  0  >-y-  h 

ges  with  the  surrounding  medium  at  temperature  T  occur  on  the  upper  r  irfa'-e 
of  the  device.  This  heat  transfer  is  considered  to  be  linear  arid  charm  ‘  ■■r'r/.it 
by  the  transfer  coefficient  M. 


2\ 


!.<•  t  0  l  x  ,  t )  be  th<*  temperature  inside  the  crystal .  The  boundary-c«,ndi t  ions  .in- 


t  !.<•  fo  1  low  i  ng 

e  (  -2h  ,  t)  -  <t>  ( t ) 

fit) 

—  +  H  (0  -  T  )  =  0  x  =  0 
dx  o 


Consequently  temperature  distribution  within  the  crystal  takes  the  farm 
from  the  Puhamel  theorem 


v  _ 0 _ 

,  (  'till  *  4hJiiJ +p,‘) 
ml  1 1 


B  B  ' 

n  ,  .  n 

?h  (X,'1,)  v?  0X1 


W  1  til 


t  f'Vt 

I  ( t )  =  /  exp  +  <)>  ( X  )  dX 

o 


(  1  ) 

(?) 

■i'T  i  7eij 

I  ;  i 

(4) 


•  is  the  thermal  diffusivity  constant  of  the  crystal  in  the  x  direction  and 

B  is  the  nth  root  of  the  transcendental  equation 
n 

B  cot  6  +  2h  H  =  0  (b) 


t’orisidering  that  <t>  ( t )  is  a  slowly  varying  function  of  time,  the  integration  of 
I  ( t  )  has  been  performed  by  parts  and  time  derivatives  higher  than  the  first  order 
are  neglected.  Then  the  temperature  distribution  is  given  by  the  relation 


i  x ,  t ) 


( B 7  + 
n 


4hJ II2  ) 


n- 1  B  (?h\)  +  4h2  H 2  +  B 2 ) 
n  n 


sin  —  (xfPh)  f  <t>  (  1  ) 


»(t) 

k  B  2 /4h J 
n 


T  if)  is  the  first  time  derivative  of  <t>(t. ). 

•  linee  the  internal  temperature  r*  ( x  ,  ♦  >  presents  a  sp.it  la  l  di.  I  r-i  to.  ’  i  ,  ih--nii.il 
■tresses  and  strains  appear. 

a"--,  mse  of  the  small  thickness  of  the  plate  with  res;  •••  t.  t  the  other  d  i  m--ns  i -uis 

1 .1 ) 

ii>  tensions  are  supported  iri  the  thickness  direction,  following  Holland  ,  the 
id  i  nbat.  i  <■  deformations  and  stresses  are  evaluated  by  taking,  into  arcom.i.  merhurw 
at  boundary  condition:;  corresponding  to  a  plate  wlucti  is  free  t.o  •‘xpainl. 


strains  and  de forma - 


Frequency  shifts  are  deduced  from  formulas  for  stresses, 
tions  by  the  method  presented  in  the  previous  1981-report.  As  in  this  calcu¬ 
lation,  summation  appearing  in  relation  (6)  is  truncated  after  the  sixth  term. 


A ui  * 

—  -  •* 

(0 

O 

■  ‘Mllpf't  .  I  t  ' i!**'  <*'!'!  i!  *  *’  '/'TO. 

If  reference  temperature  Tq  is  different  from  zero,  this  relation  is  written 
in  the  following  form: 

—  a  f  <J>  ( t )  -  T  ]  -t  a  ♦  ( t ) 
u>  o 

In  the  usual  temperature  range  from  o °C  to  10()"C  the  value  of  a  is  not  ,tc|  en- 
dont  on  the  point  of  the  frequency-temperature  characteristic,  heruu: T.r.'r; 
of  elastic  and  piezoelectric  constants  have  to  he  con:-,  i  dercd  as  constant  in 
tiie  first  order. 
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Experimental  results 


Experiments  on  the  dynamic  thermal  sensitivity  were  performed  on  different 


quartz  cuts.  SAW  delay  lines  were  submitted  to  temperature  cyclings  by  using 

„  t. _ _ _ _ i  _  (5) 


a  temperature  controlled  oven 


Fig.  2  shows  experimental  curves  obtained  by  sweeping  temperature  over  a  tempe¬ 
rature  range  of  8°C  about  10.5°C  with  sweep  periods  respectively  equal  to 
885  s  (a),  1  330  s  (b)  and  2  660  s  (c)  for  a  ST-cut  quartz  SAW  delay  line  at 


frequency  about  87.7  MHz. 


-v  -6 

The  measured  a  is  equal  to  5.6  10  s/K. 


■N, 


Fig.  2 


Experimental  sinu¬ 
soidal  temperature 
cyclings 


lernpent^r**  :r. 


wh  i  '  h  ) i.  iv<- 


Table  1  summarizes  measured  values  of  a  for  AT,  ST  arid  FST-cuts 
turnover  point  of  static  frequency-temperature  characteristics. 


1 

I  cuts 

AT 

ST 

FST 

!  a 

|  measured 

1  (s/K) 

-V  10’6 

ri 

S.b  10 

■'  in"6 

_ 

Table  I 

Values  of  the  dynamic  thermal  coefficient  a 
for  different  orientations  of  quartz  crystal 


In  temperature  ranges  beyond  the  turnover  point ,  the  presence  of  large  frequency 
shifts  coming  from  static  temperature  behavior  (non-zero  first  order  T.C.)  makes 
it  difficult  to  distinguish  between  static  and  dynamic  effects  and  thus  to  reach 
the  actual  value  of  a.  Consequently,  it  is  interesting  to  get  theoretical  values 
of  a  for  quartz-cuts  with  zero  first  order  static  temperature  coefficient. 


Theoretical  results 


The  dynamic  coefficient  a  was  calculated  for  several  cuts  of  quartz  crystal  as 
a  function  of  the  linear  transfer  coefficient  H.  To  compare  with  the  above- 
mentioned  results,  theoretical  values  of  a  for  ST  and  FST-cuts  are  reported  in 

Table  II,  the  normalized  transfer  coefficient  2hH  being  equal  to  100.  (see 
Appendix  II) 


cuts 

1  ST  ^ 

j  FST 

% 

theor. 

I  I 

*  -0  1 

-6 

<  s/K) 

-4  io  : 

1  1 

1 _  ! 

0.13  10 

i 

Table  11 

One-dimensional  theoretical  valuer,  of  a 
for  ST  and  FST  cuts 


Consequently,  the  one-dimensional  theoretical  model  with  linear  transfer  along 
the  upper  face  of  the  plate  seems  to  show  that  the  FST  cut  has  a  smaller 
sensitivity  to  dynamic  thermal  effects  than  the  ST-cut. 


hec.uise  of  the  interest  of  studying  cuts  with  zero  first  order  statu-  thermal 
coefficients,  different  double-rotation  cuts  ($,6,4»)  were  studied  and  the  t.heo- 
ret  ical  a  calculated.  The  corresponding  values  are  reported  in  Table  111.  .IT 
and  FST  quartz  are  included  for  comparison  as  is  the  bulk  SC  cut. 


Cuts 

e 

a  (s/K) 

_ 

V 

(m/s) 

Power 

Flow 

Angle 

ST 

0 

42.5 

o 

-  4.0 

io"6 

3156 

0 

FST 

f>.34 

-41.5 

26 

1.  i 

10-7 

3584 

3.8 

•l 

3 . 78 

36 

-  0.47 

Hf° 

3271 

r .  7 

77 

65.68 

-  0.65 

.o'7 

3017 

-  l.;> 

M’.W  ('(ITS 

77 

67 . 4 

27 

0 . 32 

io-ts 

3396 

-  9.8 

HO 

17 

?'/ 

-  0.6 

10"'* 

3669 

7.7 

102 

1 1  .06 

45 

-  0.12 

io-r’ 

3521 

-  1  .4 

BULK  SC 

22.5 

34.3 

3.0 

10-7 

3608 

Table  III 

Theoretical  values  of  a  for  double  rotation  cuts  with  zero  first 
order  static  thermal  coefficient  obtained  with  the  one-dimensional 
model  (angle  values  are  in  the  standard  YXwfct-notation ) 


It  should  be  noted  that  any  theoretical  value  a  is  smaller  than  1.0  10  s/K. 
Furthermore  it  is  interesting  to  point  out  that  the  cut  defined  by  angles  cf> 

=  77°,  <(>  =  65.68°  and  ^  =  36°  exhibits  a  T.C.  a  weaker  than  1.0  10-7  s/K. 
Also  it  seems  reasonable  to  search  for  cuts  with  zero  dynamic  coefficient  if 
at  the  same  time,  these  cuts  have  a  zero  first  order  static  temperature  co¬ 
efficient. 


IVVO-DIMI'.NSIONAl.  TIII.KMAI.  Dl  MI'S  ION  MODI  I 


In  order  to  model  more  completely  the  dynamic  thermal  behavior  of  SAW  devi¬ 
ces,  a  two-dimensional  thermal  diffusion  model  was  built.  In  the  axis  system 
Oxyz  the  plate  has  an  infinite  length  along  Oz,  and  a  cross  section  2i-iong 
and  2h-thick  (Fig.  3).  x 


F  i  g .  3 

Two-dimensional  temperature  distribution 


A  time  variable  temperature  $(t)  is  applied  at  y  =  ±  l.  The  boundaries,  x  *  h 
exchanges  heat  with  the  external  medium  maintained  at  T  ,  with  a  linear  heat 
transfer  such 

~+H(6-T)=0  at  x  =  +  h 
3x  O 

3  A 

-  +  H  (0  -  T  )  =  0  at  x  =  -  h 

3x  o 

The  temperature  (x,y,t)  inside  the  plate  obeys  the  diffusion  equation 

—  4.  3?e  I  i®  o  rn 

3x7  ks  3 7 y  K  3 1 

where  k7  =  X  /X  ,  X  and  X  are  the  thermal  conduct. i vif .i »•>;.:  along  1  *x  it.  I  >»v  ar.d 

x  y  x  y 

<  is  the  diffusion  constant  along  Ox. 

The  boundary  conditions  are  the  following 

tor  t  >  o  -u  v  x  <  h  y  -  *  *•  ♦<'•*  n  tJ 

V  -  l  0  <t>(  !  ) 


(7) 

(8) 
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for  t  >  0  -8.  <  y  •  8, 


x  --  -t  h 


—  ♦  H  10  -  T  )  ■-  0 
3  x  o 


x  -  -  h 


—  +  H  (f)  -  T  )  --  0 
3  x  o 


The  initial  conditions  are 

t  -  'i  f)(x,y,n)  -  1 1  -8.  *•  v  8- 

•  •  i 

-h  •'  x  '  h 

The  solution  of  equation  (9)  with  (10),  (11),  (IP)  is  calculated  in  two  steps 

a  >  first.,  by  calculating  a  peculiar  solution,  called  V(x,y,t)  obtained  with 
driving  temperature  <t>  ( t. )  considered  as  a  constant  temperature  <t> . 

b)  then,  the  solution  of  the  problem  with  a  time  dependant  excitation  <Ht)  is 
obtained  by  application  of  the  Duhamel  theorem  as  in  the  one- dimensional 
model. 

the  first  step  (a)  is  solved  along  classical  resolution  methods  of  the  diffu¬ 
sion  equation  with  boundary  conditions.  Solutions  are  developed  into  series 
such  as 
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i  is  the  nth  root  of  the  Lransceudeu  t->  I  equation 
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Solution  of  the  second  step  (b)  corresponding  to  a  variable  <t>  ( t. )  i: 


pivot,  h.y 


the  relation 
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"i’h  i relation,  Like  in  the  cade  of  the  oiie-c.  irm*n:;  iwn.i  L  model,  ha:;  two  tom;;. 
The  first  one  is  proportional  to  the  appLied  temperature  4>(t)  and  i  nt.r-  tduees 
frequency  shifts  as  in  static  thermal  behavior.  The  second  one  propoj t i oua 1  to 
the  first  derivation  of  $(t)  versus  time  is  the  term  causing  the  dynamic  ther¬ 
mal  behavior. 


CONCLUSION 


Last  experimental  results  indicate  that  SAW  devices  present  a  lower  sensitivity 
to  dynamic  temperature  effects  than  first  experiments  indicated.  In  particular 

^  _  C. 

a  dynamic  T.C.  a  for  ST-cut  equal  to  5.6  10  s/K  and  a  smaller  coefficient  for 
FST-cut  have  been  found. 

A  one-dimensional  model  taking  into  account  a  heat  transfer  function  was  used 

*v 

to  reach  theoretical  values  of  the  dynamic  TC  a  for  several  cuts  having  a  zero 
first  order  static  temperature  coefficient.  Values  of  a  are  always  smaller 

t-6 

than  1.0  10  s/K.  To  model  more  accurately  the  dynamic  thermal  behavior  of 
SAW  oscillators,  a  two-dimensional  model  of  thermal  diffusion  has  been  Pro¬ 
posed. 

In  order  to  complete  the  above  mentioned  calculations  the  next  steps  in  the 
work  concerning  a  fast  warm-up  oscillator  should  be  as  follows: 

-  calculation  of  thermal  stresses  and  strains  in  the  two-dimensional  model 

-  calculation  of  corresponding  sensitivity  in  the  case  of  cuts  with  zero 
first  order  static  T.C.  and  comparison  with  previous  results. 
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APPENDIX  I 


Experimental  temperature  controlled  oven 
for  dynamic  temperature  effect  measurements 
(fig.  4) 


For  this  type  of  measurement,  the  device  must  have  a  very  weak  thermal 
inertia.  We  have  chosen  a  temperature  controlled  oven  with  a  flat  heater. 
The  SAW  delay  line  is  held  on  a  face  of  the  heater.  Temperature  measure¬ 
ments  and  control  are  performed  with  an  H.P.  temperature  probe  which  is  in 
contact  with  the  heater  on  the  other  side. 
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